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We analyze the vector meson formulation of the BPS Skyrme model in (3+1) 
dimensions, where the term of sixth power in first derivatives characteristic for the 
original, integrable BPS Skyrme model (the topological or baryon current squared) 
is replaced by a coupling between the vector meson and the baryon current. We 
find that the model remains integrable in the sense of generalized integrability and 
almost solvable (reducible to a set of two first order ODEs) for any value of the 
baryon charge. Further, we analyze the appearance of topological solitons for two 
one-parameter families of one vacuum potentials: the old Skyrme potentials and the 
so-called BPS potentials. Depending on the value of the parameters we find several 
qualitatively different possibilities. In the massless case we have a parameter region 
with no skyrmions, a unique compact skyrmion with a discontinuous first derivative 
at the boundary (equivalently, with a source term located at the boundary, which 
screens the topological charge), and Coulomb- like localized solitons. For the massive 
vector meson, besides the no-skyrmion region and a unique C-compact soliton, we 
find exponentially as well as power-like localized skyrmions. Further, we find (for 
a specific potential) BPS solutions, i.e., skyrmions saturating a Bogomolny bound 
(both for the massless and massive vector mesons) which are unstable for higher 
values of the baryon charge. The properties of the model are finally compared with 
its baby version in (2+1) dimensions, and with the original BPS Skyrme model, 
contributing to a better understanding of the latter. 

PACS numbers: 11.30.Pb, 11.27.+d 

I. INTRODUCTION 

Among the effective field theory approaches to strong interaction physics at low energies, 
the Skyrme model [TJ plays a prominent role. The primary fields in the Skyrme model are the 
pions, whereas baryons, nucleons and nuclei are described by collective nonlinear excitations 
of the fundamental degrees of freedom of the theory, that is, topological solitons. The Skyrme 
model is very successful in the qualitative description of physical properties of nucleons and 
nuclei. First of all, the requirement of finite energy field configurations leads to an effective 
one-point compactification of the (three-dimensional) base space with the resulting topology 
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of the three-sphere S 3 . Field configurations may, therefore, be interpreted as maps from this 
base space S 3 to the field space (the group manifold SU(2)), which are characterized by an 
integer-valued topological degree or winding number. In the Skyrme model, this winding 
number is identified with baryon number, which is known to be conserved to a high precision. 
Further, a collective coordinate quantization of some light degrees of freedom (concretely, 
spin and isospin) about classical soliton solutions may be performed [2], such that baryons 
with odd baryon number are always quantized as fermions with half odd-integer spin and 
isospin, as obviously must hold true. The resulting quantum states may be identified with 
the nucleons and with both fundamental and excited states of nuclei, where the comparison 
with the experimentally measured spectra of nuclei leads to rather satisfactory results in 
those cases where a detailed calculation has already been performed (see e.g. [3]). On a 
more quantitative level, the Skyrme model, nevertheless, has some known drawbacks. First 
and foremost, higher soliton solutions correspond to rather strongly bound one-soliton bound 
states, which is in conflict with the small binding energies of physical nuclei. This problem is 
related to the fact that, although there exists a BPS bound already for the original Skyrme 
model, nontrivial soliton solutions cannot saturate this bound. The question of how to 
improve the Skyrme model towards an (almost) BPS theory is, therefore, an important 
issue. There exist two main known possibilities to improve this situation. One may modify 
the Lagrangian without altering its field contents, or one may introduce additional fields. 
The original Skyrme Lagrangian consists of two terms (the subindices refer to powers of first 
derivatives, and U is a SU(2) matrix), 

C = C 2 + £ 4 , (1) 

the so-called nonlinear sigma model term 

f 2 

£ 2 = Tr (rfdpU U ] d^U) (2) 

and the quartic Skyrme term, 

C, = -^- 2 TT([^d,U,U^ u U} 2 ), (3) 

which is required to circumvent the standard Derrick argument for the non-existence of 
static solutions. A first obvious generalization is the inclusion of a potential term 

C = -V*V(U,U*), (4) 

which is usually introduced to provide a mass term for the pions. Secondly, if a proper 
Hamiltonian formulation is required (i.e., no higher than second powers in time derivatives), 
then the only possible further generalization is provided by the following sextic term (the 
baryon number current squared) 

£ 6 = W^B", (5) 
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where is the topological (or baryon number) current 

B» = — ^-Tr (e» vprj U ] d u U U ] d p U U^d a U). (6) 

247T 2 

We remark that from the point of view of the Derrick scaling argument the sextic term is 
as good as the quartic Skyrme term. The generalized model consisting of all four terms 
has been studied and applied to the phenomenology of nucleons and nuclei, too, although 
the resulting numerical calculations are quite involved [1] - [6]. Recently it has been found 
that the restricted model consisting only of the potential and sextic terms (the so-called 
BPS Skyrme model) has a BPS bound and exact soliton solution saturating this bound [7J, 
[8]. These classical solutions, therefore, correspond to nuclei without binding energies and 
realistic, small binding energies may be introduced both by quantum corrections and by 
small contributions of additional terms [5J- The BPS Skyrme model has further intriguing 
mathematical properties like, e.g., infinitely many symmetries and conservation laws as a 
consequence of its generalized integrability [TO]. These symmetries contain the symmetries 
of an incompressible ideal liquid and allow, therefore, to reproduce some features of the 
liquid drop model of nuclear matter. 

Another possibility to overcome the shortcomings of the Skyrme model consists in the 
inclusion of further (e.g., vector) fields. Firstly, one may couple the Skyrme fields to the 
electromagnetic field. The correct electromagnetic coupling of the Skyrme model was first 
derived in and soliton solutions of the resulting Maxwell-Skyrme system have been 
studied in [12], [13] (the equivalent problem for the baby Skyrme model in one dimension 
lower has been studied in [H], [15], and, for the BPS baby Skyrme model, in [16J). Secondly, 
one may couple the Skyrme fields to vector mesons. In a recent investigation, a Skyrme 
theory coupled to an infinite tower of vector mesons was derived from an instanton holonomy 
in one dimension higher, where the exact BPS property of the resulting Skyrme theory may 
be traced back to the self-duality of the instantons [17] . The infinite tower of Kaluza-Klein 
vector modes induces a flow to a conformal BPS theory which is a Minkowski space version 
of the Sakai-Sugimoto model (for further developments and applications to baryonic physics 
see, e.g., [H]). So one main difference between the BPS Skyrme proposals of [T7J and of 
[7] resides in the different symmetries which are associated to the BPS property (conformal 
symmetry in the former case, volume-preserving diffeomorphisms in the latter case). Also, 
it is not obvious how to include a potential term into the conformal setting of [17], and this 
issue has not yet been completely resolved, to the best of our knowledge. Any truncation to 
a finite number of vector mesons of the theory of [IT] leads to a theory which is no longer 
exactly BPS, but has rather small binding energies. 

Among all the vector meson couplings to the Skyrme model, there is one which results in 
a Lagrangian which is rather similar to the sextic term (the baryon current squared) above, 
namely the so-called omega meson u^. Indeed, the omega meson couples to the topological 
current, C ~ ^-B^, and upon integrating out the omega meson, the sextic term is recovered 
in the limit of infinite vector meson mass. The importance of the omega meson is also related 
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to the fact that it transfers the physical effects of the chiral anomaly to baryonic matter and 
that it prevents the solitons from shrinking due to the appearance of a short-range repulsion 
in nuclear interactions [T9]. Moreover, it blocks the flow to a conformal theory. 

We remark that the omega meson may be integrated out also for finite or zero meson 
mass, resulting in a nonlocal effective self-interaction of the topological current with an 
integral kernel K^ u of the Yukawa or Coulomb type, 



although we shall not pursue this approach in the present paper (i.e., we will always main- 
tain the omega meson explicitly). The same nonlocal interaction induced by a Yukawa or 
Coulomb integral kernel has been investigated for the nonlinear Schroedinger equation in 
lower dimensions in [20] . The scaling behavior of the vector meson terms (the coupling term 
and the standard kinetic term) is, in fact, such that it stabilizes the solitons without the 
need for a Skyrme term [21] - [23] , and the resulting theory consisting of £ an d C-2 coupled 
to the omega meson has been studied recently in [21], where it was found that the soliton 
solutions are quite similar to the solitons of the Skyrme model. 

It is the purpose of the present paper to study in detail the vector meson version of the 
BPS Skyrme model, which is obtained from the vector model described in the preceding 
paragraph by suppressing the sigma model term £ 2 , and to compare its properties to the 
standard BPS Skyrme model. First of all, we shall find that the infinitely many symmetries 
of the BPS Skyrme model carry over almost unaltered to the vector BPS Skyrme model 
and, consequently, we will be able to perform most of the calculations analytically, like in 
the BPS Skyrme model case. Concerning soliton solutions, we shall find that the solitons 
of the vector BPS Skyrme model are rather different from the ones of the standard BPS 
Skyrme model, at variance with the results of [21] for the full Skyrme model and its vector 
version. This different result is, in some sense, expected, at least for certain potentials. The 
reason is that for potentials which include a pion mass term, the suppression of the term £2 
corresponds to the limit of infinite pion mass in the sense that linear fluctuations of the pion 
field are completely suppressed. On the other hand, physically the pions are the lightest 
effective particles. Hence, the inclusion of vector mesons with a finite mass and a standard 
kinetic term in the BPS Skyrme model in some sense reverses the typical mass hierarchy 
of low-energy QCD. It would, therefore, be surprising and cast some doubt on the viability 
of the BPS Skyrme model as an effective theory for strong interaction physics, if it lead to 
qualitatively similar solitons as in the vector model with its inverted mass hierarchy. We 
remark that the vector versions of the baby Skyrme model and its BPS restriction have been 
investigated in [25J and in [26J, respectively. These results are extended and generalized to 
the 3+1 dimensional situation in the present paper (for the BPS case), and in [24J (for the 
full Skyrme model case). 
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II. THE ui- VECTOR MODEL 

The vector version of the BPS Skyrme model is given by the following Lagrange density 
£ = -fj?V(U, tf) - ^{d^ - d u u,f + l -M 2 ul + X'upB" (7) 

where is the baryon current. 

We use the standard parametrization of the SU(2) chiral field 

U = e* a ' 9 = co&€ + ifan£n-a, n 2 = 1, 

where f are the Pauli matrices, £ is a real field and n is a unit three component vector field, 
which is further related to a complex field u by means of the stereographic projection 

ft = ; — ttt (u + u, —i(u — u), 1 — \u\ 2 ) . 

1 + \u\ 2 ' 

Then, 

C = -» 2 V(U, C/t) - \{d^ v - d u u,) 2 + \M 2 ul + ^^ u^%u a u (8) 

where A is a new constant related to A'. The potential term is assumed to depend only 
on TrU i.e., on the scalar field £. Concretely we are going to analyze in detail a family of 
potentials which provides a generalization of the usual Skyrme potential 

V= (^^) Q = (l-cosO°. (9) 

Soliton solutions of the Skyrme model with the old potential (a = 1) have been studied, 
e.g., in [27], [28], whereas generalized potentials were investigated, e.g., in [29], [30]. The 
pertinent field equations take the form 

d,F^ + MV + J A + S ^5 )2 6^£ a ^n 7 = (10) 

iXe^d^pu-, = (11) 
+ = 0- (12) 

[i + \u\ J 

We assume the natural static ansatz 

u = co = co(r), £ = £(r), u = v(6)e m<t> (13) 
and the other u-i = 0. Then, their static versions are 

V 2 r u - M 2 u = iX , 2 , 2 V r £(V^V M - V e uV^u) (14) 

[l + \u\ ) 
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iX sin £ 
1 + Id 2 ) 2 



VruiVeuV^u - VeuVtpu) + = 0. (15) 



Observe that one of the equations of motion i.e., (11), is obeyed identically by the ansatz 
without any restrictions on the form of the ansatz functions. However, the ansatz is compat- 
ible with the remaining two equations only for a very restricted form of the complex field. 
Namely, we need that 

(TTR^ (6) 

is a function of r only and does not depend on the angular variables. Then, the remaining 
equations (14), (15) become ODEs depending entirely on r. So, for our ansatz 

VeuV^u — VeuV^u _ 2in Wg 
(1 + Id 2 ) 2 = ~r 2 sin#(l + i; 2 ) 2 



(17) 



where the function v must obey the appropriate boundary conditions such that u covers 
the target space two-sphere at least once. The well-known solution which covers the target 
space latitude exactly once is 



9 VeuV^u -VeuV^u in 

v (0 = tan - =>- -r- — , , nso — — = -. 18 

2 1 + M 2 2 2r 2 



Inserting this result into the static equations we get 



V 2 w - M 2 u = n\— ^V r £ (19) 



2r 2 



nX sin 2 £ 



V r uJ + yrV^ = 0. (20) 



2r 2 

Further, we may use the last formula to eliminate the derivatives of the vector fields from 
the first equation. Then we arrive at a set of two first order ODEs of the gradient flow type 

1 a ( 4 2fi 2 V i nA, 



r 2 \ nX sm £ 4 ' 



Ur^-r^JL^L-. (22) 
nX sm 2 £ 



The fact that the static equations of motion can be reduced to a set of two first order ODEs 
(solvability) may probably be related to the existence of infinitely many conserved charges 
(integrability) . This set of equations must be equipped with the proper boundary conditions 
which guarantee nontrivial topology, 

£(r = 0) = 7T, £(r = Rq) = 0, (23) 

where Rq can be finite (compactons) or infinite for usual solitons. Further, the boundary 
conditions for the vector meson field are 



Ur {r = 0) = 0, u(r = R ) = 0. 



(24) 
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Notice that the last formula leads to the following condition for the behavior of the profile 
function at the origin 

*™ (^^TT-r) = 0. (25) 
sin 4l r J / 



r->0 



A. Massless case 



For the massless vector meson field one may integrate the first equation of motion com- 
pletely 

,9„ 2 y, 

+ — (£- sine cos £) = C. (26) 



4 2u 2 Vc nX 

r 2 — _|_ ' 

nX sin 2 ^ 4 



The integration constant C can be easily determined using the assumed boundary conditions 



and (25) 



£7==* (27) 

Thus, for the massless case we have obtained exact solutions (although usually it is not 
possible to write them in a closed form, i.e., as £ = £( r )) for any value of n. The questions 
whether these configurations correspond to a nontrivial topology and how the solutions are 
localized, are determined by the particular form of the potential. 



1. No solutions - a € [1, 1) 

For the most interesting case of the old potential the profile equation reads 

4 2u 2 1 n\.^ . ^ ^ nnX , . 

i JJ Y—^ + "T~(£ - sine cos = -p. 28 
nX sine 4 4 

However, the profile function defined by this equation cannot reach the vacuum e = for 
any value of r ^ 0. Indeed, the left hand side is singular at such a point while the right 
hand side is obviously finite. Thus, there are no topologically nontrivial configurations for 
the old potential. This can be generalized to all potentials with a < |. Then, the left hand 
side of the profile equation still is singular at the vacuum e = 0, 

cos | V 2 / 4 4 



In this case we get 



2. Compacton - a = | 



4 ^2/i 2 1 nX 7mA 

r — j + —(e-smecose) = — r- 30 

nX cos | 4 4 



8 



which at the vacuum value of the profile function takes the form 

4 yV _ ™A 

°^a~ " T- (31) 

Obviously, we find a compact skyrmion for which the vacuum value must be reached at the 
finite radius 

As in the case of the vector BPS baby Skyrme model [26], any compacton solution is, in 
fact, a solution in the space C of continuous functions, but not in the space C 1 of continuous 
functions with a continuous first derivative, because the derivative is discontinuous at the 
compacton boundary. Equivalently, any compacton is a solution of the field equations with a 
Dirac delta source located at the boundary of the compacton, which screens the topological 
charge generated in the interior of the skyrmion, 

-«9 r (rV + — (£-sm£cosOJ = —5{R-r). (33) 

Hence, the total charge of this configuration is zero. 



3. Coulomb-like localized solutions - a > | 



Now we find a power-like tail for the profile function of the skyrmion, where the power 
depends only on the potential but not on the topological charge 



7in 2 X 2 



2a-Z I 



f 2a-3 

while the vector meson field is universally of the Coulomb type 

7mA 



0J r 



4r 2 



+ 



The corresponding energy reads 



E = 4tt / r 2 dr 
Jo 

POO 

4tt / r 2 dr 
Jo 



An 
2T 2 



uj sin 2 ££ r 



"°° d 

— 47r / — (r 2 uu r ) 



o 



dr 



An I r 2 dr 
'o 



"V(0 + -oj 2 . 



(34) 



(35) 

(36) 
(37) 
(38) 



where, using the equations of motion, we have combined the second and third term into 
a total derivative plus an additional term. Since, due to the boundary condition and the 
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asymptotical form of the solution, the total derivative term is zero, we arrive at a posi- 
tive definite expression for the total energy. Obviously, the second term is localized like a 
Coulomb electric field. On the other hand, the potential term converges in a way which 
depends on a. Indeed, for £ — > 

rV(0| h „ = r 2 (l - cosOl „ „ ~ r 2 e a \ h „ ~ r 2 ' ' A ' / ' A " 

y ^'\onshelL * \onshell ' \onshell 

It approaches the Coulomb like localization for a — >• oo while for any other a > 3/2 we get 
a stronger convergence. 



B. Massive case 



In the massive case, the field equation for the vector meson can be rewritten as 



<9, 



-(£ - sin£cos£) 



-M z u 



nX sin £ 4 
where, for simplicity, we introduced a new variable x = r 3 /3. Then, knowing that 

2{j?a (1 - cos£) a_1 



UJ, 



sin£ 



and acting with d x on (|39|) we arrive at 
(3x 



0: 



sin£ 



n 2 A 2 



8/i 2 a 



(f - sin £ cos f ) 



M 



^l-cosQ"- 1 
sin£ 



(39) 



(40) 



(41) 



In order to classify skyrmion solutions from the point of view of their asymptotic behavior 
it is enough to analyze this equation in the vicinity of the vacuum value £ = 0. Then we get 



(3x) 4 / 3 £ 2Q - 3 + 



n 2 A 2 2*- 3 
[i 2 a 3 



MX 



2 tf 2a-3 



(42) 



Performing the series expansion we find three types of possible solutions: a compacton, and 
exponentially as well as power-like localized skyrmions which are analyzed in the further 
part of this section. 

Finally, let us notice that using the field equations the total energy may be rewritten in the 
following form 



E = An r 2 dr 
Jo 

POO 

= An I r 2 dr 
Jo 



fj? (1 - cos £) a + -lu 2 + -MW 



'fi 2 (2h) a + \u 2 r + l -MW 



where 



h(r) 



1 - cos£) 



(43) 
(44) 
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1. No solutions - a < | 

For £ close to the vacuum we get 

d 2 x [(3x)^ 3 e a - 3 ] = M 2 e a ~ 3 (45) 

where both the left and the right hand side diverge. Then, the asymptotic solution (x — > oo) 
is 

which means that £ is exponentially localized. However, it leads to an exponentially divergent 



vector meson (40) and therefore to infinite energy. Hence, there is no skyrmion in this 
parameter range. 



2. Compacton - a = | 

In this case one can find a compact solution 

f(x) x<X 

M { (47) 
x > X 

tu(x) x < X 

,cU) { (48) 
x > X 

where X is the (third power of the) radius of the compacton at which / and u i.e., the 



solutions of (|39j), (40), are joined with the vacuum configuration. However, u x has a jump 
at the boundary as to x {X) = — 3/j 2 /n\\/2. Hence, the derivative of it produces a delta 
like term located at the boundary. More precisely this boundary source term cancels the 
topological charge inside the compacton leading to a zero total charge configuration, exactly 
as in the massless case. 



3. Exponentially localized skyrmions - a € (f,3] 



Assuming a < 3 one gets that £ > £ for small £. Then, the first term in (42) becomes 
the leading term under the additional condition (which must be verified at the end) that £ 
goes to sufficiently fast. Hence, we find 

d 2 x [(3x) 4 / 3 £ 2Q - 3 ] = M 2 £ 2Q - 3 (49) 

which is a linear equation for £ 2a ~ 3 . It possesses an exponential solution 




O 1 r 

1 2 3 4 5 6 

FIG. 1: Solutions with their derivatives and the energy density for the a = 2 massive case with 
the constants f/, = M = \ = n = l. 



Further, also the meson field and the energy density are exponentially localized. 

An example of a solution of this type has been calculated numerically for a = 2. Then, 

performing the expansion around the center for the equations, we find (we use again the 
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field variable h instead of £) 

h ir)^ + ^- M2 ^- 5 ^^ ^ + ... (51) 



2/i 2 y^ 2 ^ 2 M 2 /i 2 + 6n 2 A 2 /f 4 

w(r) ~ i; + — ^ — r - ^ x ^ 9 ; „ 2 ^ 2f2 r + ••• (52) 



n\f 2 lOnAv^/a /i 2 + n 2 X 2 f, 



where f'2 and Vq are the free parameters we can vary in order to find a numerical skyrmion 
solution with the correct asymptotic behavior for large x (exponentially decreasing) via a 
shooting from the center. For A = jj, — M — 1 and the simplest charge one skyrmion they 
take the following values 

f 2 = -1.8384364, v = 1.14, (53) 

and we find the solutions presented in Figure [Tj The static energy corresponding to the 
energy density of Figure [T] is E = 5.73645. 



4- Power-like localized skyrmions - a > 3 



For a > 3, the second term in (42) is dominating which leads to a different type of 
vacuum approach. Now we have 

2x2 2 "-3 

71 A - d 2 e = £ M (54) 

which gives a power-like localization 

1 

1 \ Q ~ 3 



■ ( 55 ) 



This behavior is reproduced by numerics. For example, for a = 4 after expanding around 
the center we find 

h(\ i-u* 2 8M 2 /. 2 - 200/i 2 / 2 + n 2 A 2 /| 4 

W ~ 1 + hT S^ + r^A 2 / 2 ) hT + (56) 

u(r) ~ „ - 16 ^ ^ _ 8 8MV + 13n 2 A V/| 
The free parameters are now fixed as 

/ 2 = -4.98624, w = 1.64, (58) 

which leads to the solution presented in Figure [2] with the static energy E = 7.53696. 




1 2 3 4 5 6 

FIG. 2: Solutions with derivatives and the energy density for the a = 4 massive case with the 
constants fi = M = X = n=l. 



C. BPS case 



Interestingly, there is a potential for which the energy density can be written as a to- 
tal derivative. This potential possesses only one vacuum, where it looks similarly to the 
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previously considered case a = 4. Globally, however, it has the following form 



Vbps = ^ (£- cos 4~ sin f) 2 
It has the nice property that it allows to simplify one of the field equations 



nA 



sin 2 £u r = —fi 2 (£ — cos £ sin £) sin 2 £ 



nA 



-pt 2 (£ - cos£sin£) 



or, using the variable x = r 3 /3 introduced in Section II. B, 

nX 2 . . 

—u x = -fi (4 - cos 4 sm£) . 

Now, using the equation for the meson field, one can compute the total energy 



E = 4tt 



r <ir 



-IV 



47T 



47T 



47T 



<ix 



dx 



An 
2r 
An 
4r 
An 



2 o;sin ££ r 



At Vbp S (£) - — wsin 



4 



nA 
V 



n 2 A 2 
16> 2 



w 2 + 



tUJU 7 



where the last step follows from the derivative of (61), i.e., 

—u xx = -2/i sm 44V 
Then, combining the two terms into a total derivative we arrive at 



2 \ 2 r°° 



E = 7T- 



n 2 X 



dx [ul + uu xx ] 



2 \2 r°° 



7T- 



n^A 



n 2 A 2 



<ix 



4/i 



(59) 
(60) 

(61) 

(62) 
(63) 
(64) 

(65) 
(66) 



dx—(uju x ) = -iv——z-u)(0)u) x (0) (67) 



where the asymptotical vanishing of the meson field has been used. Further, the value of oj x 
at the origin is fixed by equation (61) and the boundary condition for the profile function 

m = * 

27T/Z 2 



w*(0) = — 



nA 



Hence, finally 



E 



ir 2 nX 



w(0). 



(68) 



(69) 



Observe that this formula is valid for both the massless and the massive case, because the 
vector meson field is sufficiently localized also in its massless version. In the massless case, 
it is possible to find the value of the meson field at the origin analytically. The massless u 
meson equation (26) has the solution 

4" — cos £ sin £ 



7T 



1 + ^r 4 ' 



(70) 



-i 2 A 2 ' 
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Thus, 

nX 7ifj, 2 r 2 



* w ' = "rrjp (71) 



which may be integrated to the exact expression 



3 /2 

x V2 UVfyJ 



27r/i 2 1 / nX \ 



nX 



In -= I — 2 arctanfl — y/2ar) + 2 arctan(l + y/2ar) — 2n 

\ar 2 + y/2ar + 1 J 

where a = 2y/2fi/nX. Hence, its value at the origin is 



(72) 



/ x 7T 2 / nX \ , . 
= T • ( 73 ) 



4 V2v / W 
Finally, the energy 

£ M=0 = ^ — Aw — ■ n5 (74) 
8V2V2 V ^ 

Obviously, a sum of n charge one solutions is lighter than the charge n soliton. Therefore, 
higher charge skyrmions are unstable although they saturate the BPS bound. Whether this 
is a consequence of the assumed spherical symmetry or a general property of this BPS model 
remains to be checked. 

In order to (almost) analytically find the value of the meson field at the origin for the 
massive case we use its field equation 



^((r 4 ^ + ^)(e-8ineco80) =-MV (75) 

Then, we switch to the coordinate x and differentiate the resulting equation by d x . Using 
(61) we arrive at 

^(((3x) 4 / 3 ^ + ^(e-sinecose)) =^^(£-sin£cos£) ( 76 ) 
which is a second order but linear differential equation for h — £ — sin£ cos£ 

d °({ {3x)4/3+ w) h{x) ) =M2h{x) (77) 

The value of the vector meson field at the origin may be related to the derivative of h by 
using formula ( 75 ) 

n A 

"(°) = -4^(0). (78) 

Unfortunately, the equation for h can be solved only numerically, despite its linear character. 
We solved it by shooting from the center. We assume the standard values for the constants 
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/j, = M = X = 1. The range for the numerical solutions has been chosen as R = 100. 
Moreover, it has been checked that the BPS energy (computed for each skyrmion with 
Q G [1, 10]) gives exactly the same value like the one we get from the usual expression 



E = Air dx 



o 



(79) 



It is clearly visible in Figure [3] that the higher charge skyrmions are unstable although they 
saturate the Bogomolny bound. Again, we observe a faster than linear growth of the energy 
with the topological charge. In fact, one can fit a curve E = a ■ \Q\ b and gets a = 5.13 ± 0.07 
and b = 1.196 ± 0.007. Although we do not have any analytical argument, the power 
parameter b is equal (within the error) to |. In any case, this dependence is slightly weaker 
than in the massless case. An example of solutions with Q — 1 is presented in Figure 
|4j Taking the large x limit, one can find that the solution possesses an exponential tail. 
Configurations for other topological charges look very similar. 

E/n 

8 - . 



4 - 



2 - 



n 

10 



FIG. 3: Total energy over the topological charge as function of n. 

Let us remark that for the family of potentials 

V = {VbpsY (80) 



built on the BPS potential (59), we observe the same pattern of solutions as in the previously 
analyzed case (the family of powers of the old potential). In the massless version, there are 
three possibilities: no skyrmions for (3 < |, a compacton with a screening source at the 
boundary for /3 = |, and Coulomb- type solutions for /3 > \. In the massive version we 
found four cases: no skyrmions for (3 < ^, a compacton (which again requires a screening 
source) for (3 = |, and two infinitely extended types of solutions - exponentially (f3 G (|, 1]) 
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h,h' 

3 - 




FIG. 4: Solutions and energy density for the BPS potential with the constants fi = M = A = 1, 
and winding number n = 1. 



or power-like localized (/3 > 1). In fact, this behavior repeats for any one- vacuum potential 
with the power-like approach to the vacuum at £ = 



V ~C 
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Here, a = 2a for the old Skyrme potentials family or a = 6/3 for the BPS family. 

All these results should be compared with the qualitative properties of skyrmions in the 
BPS Skyrme model for the potentials considered above. This is analyzed in the next section. 

III. SKYRMIONS IN THE BPS SKYRME MODEL 

The BPS Skyrme model is [7] 

C BPS = \ 2 n 4 Bl-^V(U,U^), (81) 
which, using the field decomposition introduced previously, reads 

£bps = - (1+^)4 (^6^) 2 - » 2 m- (82) 

Assuming exactly the same Ansatz for the Skyrme field, we arrive at a differential equation 
for the profile function £ 

n 2 A 2 sin 2 £ / sin 2 £ £ r 



d r ) ~ M = 0. (83) 

This equation can be simplified by introducing the new variable z (up to a numerical factor 
it is our previously defined x) 



V2fi 



?> 



Then, 



and it may be integrated to 



3 \n A 



sin 2 e^(sin 2 ee 2 ) - ^ = 0, (85) 
\sm^e z = V(0, (86) 



A. The old Skyrme potentials 

For the old Skyrme potentials it gives 

±-^sin 2 a 2 = (i-cosO Q/2 . (87) 

In the generic situation, this equation is solved by a combination of hypergeometric functions 
which is not quite illuminating. However, some general observations can be easily made. 
First of all, performing an expansion in the vicinity of the vacuum value one can identify for 
which potential compact skyrmions occur. Namely, if we assume that £ ~ 0, then the BPS 
equation at the leading order reads 

l-a 



±£ 2 = 2^r _2 (88) 
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which possesses the obvious solution 

(z-z )^ a e [1,3) 
f ~ { e--2 z a = 3 (89) 

(I)- 3 «>3. 

Thus, for a < 3 the resulting skyrmions are of the compact type. This means that the 
vacuum value is reached at a finite distance, the radius of the compacton. It is in agreement 
with the result for the old Skyrme potential (i.e., a = 1), for which compact configurations 
have been found previously. For a = 3, we have an usual exponentially localized solution, 
while for a > 3 a power-like approach to the vacuum is observed. One can find exact 
solutions which confirm this result. 

a = 2 The pertinent solution reads 

e = o z>f 2 . 

As expected, we get a compact configuration. This solution has a jump of derivatives 
at z = 7r/v^2- In this case, the jump is finite, and, as in the case of the standard Skyrme 
potential, this jump is immaterial for physical quantities like the energy density or the 
topological charge density. 

a = 3 The solution is given by the following expression 

cos - + In tan - = . 

2 4 2 

Notice that this profile function is not of compact nature but is non-zero for all z. The 
vacuum value is approached asymptotically at infinity. Moreover, this solution and its 
energy density are exponentially localized. 

a — 4 Now, the solution fulfilling the assumed topological boundary condition is 



{ + 2coti = ^( 2 + ^). 



This solution is non-zero for all z and i' z {z) — > 0, if z — > oo. It is localized like a 
polynomial in inverse powers of r. 

a = 6 For this value we derive again a solution which is non-zero for all z, 

£(z) = 2 arc cot \J$y/2z. 

It easy to see that this configuration as well as the corresponding energy density are 
polynomially localized in inverse powers of r, as for the a = 4 solution. Indeed, 

e(z) = 16V2nn\n — — 5-. 

l + (3V2z)l 
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Finally, we can compute the total energy of these solitonic configurations for any a value 
using the BPS property of the solutions. Then, in the z variable 

E = ATc\n\\fi J dzt; z sm 2 £W = An\n\\fi d£sin 2 £(l - cos£)t (90) 

= 47r 3/2 2 l +f r (l + f) |n|A (Q1) 

r (3 + f ) 

The existence of compact, exponential and power-like skyrmions in the BPS Skyrme models 
is an expected phenomenon as it has its counterpart for a lower-dimensional version of the 
model - the so-called BPS baby Skyrme model. 



B. The family of BPS potentials 

Now we consider the family of potentials constructed from the previously introduced BPS 
potential 

V = {V BPS f = Q(£ - cos£sin0 2 ) . (92) 

Locally near the vacuum, the family of BPS potentials looks exactly as the family of the 
old Skyrme potentials with the identification a = 3(3. So, the approach to the vacuum is 
the same. On the other hand, the behavior at the origin where £ = n is different. However, 
this fact only weakly (quantitatively) influences the solutions. Hence, we may consider this 
family of potentials as an approximation of the more standard ones which allows for exact 
solutions for all values of the parameter /3. 
Now the BPS equation reads 

±-J=sin 2 ^= Q(e -cos^ sin 0^ (93) 

It can be further simplified if we define a new target space variable 

77 = ^(£-cos£sin£) (94) 

which interpolates between 77(0) = | and 0. Observe that using this variable V = r/ 213 , and 
the profile equation takes a very simple form 

± ^= Vz = rf. (95) 

It can be easily solved leading again to three types of skyrmions. 
For j3 G (0, 1) we find compact skyrmions 

v= ! {V2(l - (3)(z - z))^ ze[0,z ] (Q6) 

(0 z > Zq. 
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For p — 1 we have a unique exponentially localized solution 



V = |e-^. (97) 



For /3 > 1 we find power-like localized skyrmions 



V ( .v / 2(/3-l)(^ + ^o),' ' ^ 



Here 

' = ~y/2\0-l 



Zo ~ -4-11 (I) 



Finally, the corresponding energy is 



£ = Att 1 -^-^ - . (99) 



1 + /3 V2 

IV. INTEGRABILITY AND CONSERVATION LAWS 

As one might expect, the vector BPS Skyrme model is integrable in the sense of general- 
ized integrability. In particular, there is an infinite family of conserved currents 

* = fu ^ - fu^ (100) 

where 

W = e a ^u a ^u u , ^ = e a ^u a ^u u (101) 
and G = G(u,u,£) is an arbitrary function of the target space coordinates. Then, 

d»J° = G m u^ + G uu u^ + G u d^ - G ua u^ - G uu u^ - G u d^ (102) 

+ G^£" - G^/C" = (103) 
where one has to use identities obeyed by /C M 

VC* 4 = ^/C" = 0, u l X» = u l X» (104) 

as well as the field equation 

fy/C" = 0. (105) 

Observe that this is exactly the same family of conserved quantities as for the original 
BPS Skyrme model. The interaction with the vector mesons does not spoil the generalized 
integrability property, which may perhaps be responsible for the solvability of our model. 
The set of conserved currents is even bigger if one considers the massless version of the 
model. Now, we can construct the additional family of currents 

jff = H(uu)F IM ,(uu v + uu v ) = H{uu)F^d u {uu) (106) 
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where H depends now on the modulus of the complex scalar. Then, 

d"jn = HF^d u (uu) + H'F^{uu)d u {uu) + H(d»F^)d u {uu). (107) 

The first two terms vanish due to the contraction of the antisymmetric tensor with two 
symmetric ones. Using the field equation for the mesons we get 

&>tf = -iXHe^ ^ j + uti") = 0. (108) 

This family of conserved currents is identical to the one found in the baby version of the 
model [2"o] . 

V. SUMMARY AND CONCLUSIONS 

In the present paper we investigated the vector BPS Skyrme model and its soliton solu- 
tions. In this model, the sextic term of the BPS Skyrme model is replaced by a coupling of 
the baryon current to a (massive or massless) vector field (the omega meson). First of all, we 
found that the resulting model is still integrable in the sense of generalized integrability [10] 
and possesses infinitely many symmetries and conservation laws, which are, in fact, identical 
to the ones of the original BPS Skyrme model. Further, we found that, for a spherically 
symmetric ansatz, the static field equations may be brought to a first order form of the 
evolutionary (or gradient flow) type, where this fact is probably related to the integrability 
properties of the model. For a massless vector meson, the resulting equations can even be 
integrated completely, although the solution can usually be given only in an implicit form. 

The soliton solutions of the vector BPS Skyrme model are quite different from the corre- 
sponding solutions of the original BPS Skyrme model. For the standard pion mass potential, 
for which the BPS Skyrme model has a compacton solution, the vector model has no soliton 
solution at all, because any local solution cannot be extended to a global one with the correct 
boundary conditions for a topological soliton (i.e., all formal solutions have infinite energy). 
As explained already in the introduction, this difference has to be expected, because the 
vector BPS Skyrme model with a pion mass term corresponds to a situation where the mass 
hierarchy is reversed with respect to the case of low energy QCD. This seems to imply that 
if one wants to include omega mesons in a physically reasonable way into the BPS Skyrme 
model (i.e., without altering the physical mass hierarchy), then the limiting case of infinite 
meson mass should be considered. But this brings us back exactly to the original BPS 
Skyrme model, up to a redefinition of its coupling constants. The BPS Skyrme model itself 
may, therefore, be interpreted as a vector meson Skyrme model in a certain limit, which 
might explain its BPS property from a slightly different perspective. 

For other types of potentials, with a faster than quadratic approach to the vacuum, 
soliton solutions of the vector model may exist, but their behavior is quite different from the 
normal BPS Skyrme model case. Genuine compacton solutions, which are quite typical for 
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the BPS Skyrme model, do not exist in the vector version. They only exist (for potentials 
with a cubic approach to the vacuum) for a modified field equation where an inhomogeneous 
delta function source term effectively screening the topological charge is introduced at the 
compacton boundary In other words, the vector model has compacton solutions only in 
the space C of continuous functions, but not in the space C 1 of continuous functions with 
a continuous first derivative. For potentials with an even faster (than cubic) approach to 
the vacuum, solitons with a power-like or exponential decay at spatial infinity can be found. 
Concretely, for a massless vector meson we find that the vector meson term in the energy 
density always localizes like a Coulomb term, whereas the (power-like) localization properties 
of the Skyrme field depend on the potential. In the massive vector meson case, depending 
on the vacuum approach of the potential, both terms may decay either exponentially or in 
a power-like fashion, as summarized in Fig. [5] below. 



BPS Skyrme 

a=l a=3 



C'-compacton 

exponential power-like 



Massless vector 

a=l a=3/2 a=3 




no solutions 

C-compacton Coulomb-like 

Massive vector 

a=l a=3/2 a=3 



4 f>,V i 



no solutions exponential power-like 

C-compacton 



FIG. 5: Comparison of types of solutions in the BPS baby Skyrme mode (the old skyrme potentials) 
with its massless and massive vector counterparts. For the family of the BPS potentials the picture 
is exactly the same provided a substitution a — > 3/3 has been made. 



Next, we considered the case of a specific potential (the BPS potential of Section II. C), 
where the energy density for spherically symmetric solutions may be written as a total 



24 



derivative, such that these solutions saturate a BPS bound. It turns out that the energy 
of the BPS bound grows faster than linear in the topological charge and, therefore, higher 
charge solitons are unstable w.r.t. decay into smaller ones. At the moment it is not known 
whether this is just a property of the sector of spherically symmetric solitons, or whether this 
instability (and maybe even the BPS bound) continues to hold for the full model. Probably 
a full three-dimensional numerical simulation would be necessary to clarify this issue. In 
any case, this behavior is, again, completely different from the case of the BPS Skyrme 
model, where already the spherically symmetric solitons saturate a BPS bound linear in the 
topological charge. This linear behavior is, in fact, quite important for its applications to 
low-energy QCD and nuclear physics. For further recent results on nontrivial BPS bounds 
we refer, e.g., to [21], [15], [25] • 

If we compare our findings with the analogous results for the vector BPS baby Skyrme 
model in one dimension lower [25], then most of our results are quite similar. One minor 
difference is that in the baby case, the compacton solution with the effectively screened 
topological charge appears already for potentials with a quadratic approach to the vacuum 
(i.e., with a mass term). In both cases, genuine compactons do not exist for the vector 
models. This absence of compactons is probably related to the presence of a Coulombic (or 
Yukawa-like) term in the action of the vector models, which implies long-range interactions. 
Only if the topological charge (which plays a role analogous to the electric charge due to its 
coupling to the vector meson) is effectively zero because it is screened by a source term, these 
Coulomb or Yukawa long range interactions may be absent. It is interesting to observe that 
when the electromagnetic field is coupled to the baby Skyrme field in the standard, minimal 
way, then the resulting gauged BPS baby Skyrme model allows for compacton solutions [15] 
without problems. In this case, however, the electric field (and, therefore, the electric charge 
density) must be zero for finite energy configurations, and only magnetic fields are allowed. 
There is, therefore, no Coulomb term which would give rise to a long-range interaction. 
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